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GROUND  MOTION  EFFECTS 
FROM  NUCLEAR  EXPLOSIONS:  A  REVIEW  OF 
DAMAGE  EXPERIENCE  AND  PREDICTION  METHODS 


Abstract 

# 


Ground  motion  from  nuclear  explosions 
can  damage  one-  and  two-story  buildings 
when  surface  velocities  are  as  low  as  0.1 
to  0.2  cm/sec,  while  high  explosive 
blasting  experience  had  previously  sug¬ 
gested  a  threshold  of  between  5  and 
10  cm/sec.  Subsequent  studies  indicate 
that  this  damage  correlates  best  with  the 
average  amplitude  of  the  pseudo  absolute 
acceleration  (PSAA)  in  the  0.05  to  0.2  sec 
period  range.  Present  data  indicate  a 
damage  threshold  of  about  10  cm/sec^, 
while  at  PSAA's  of  about  1000  cm/sec^, 


damage  to  over  half  of  the  buildings  can 
be  expected.  This  damage  is  architectural 
in  nature,  consisting  predominantly  of 
cracked  interior  plaster,  cracked  brittle 
masonry  walls,  and  damaged  chimneys; 
about  two-thirds  of  the  Rulison  damage 
was  of  this  nature.  The  use  of  this 
damage  prediction  method  requires  the 
correct  prediction  of  response  spectra. 

The  accuracy  of  such  predictions  has 
been  greatly  improved  by  employing 
frequency  dependent  yield  and  depth 
scaling. 


Introduction 


From  the  inception  of  the  Plowshare 
Program  it  has  been  recognized  that 
understanding  and  predicting  the  ground 
motion  generated  by  the  explosion  and  any 
associated  building  or  structural  damage 
was  necessary  if  nuclear  explosions  for 
peaceful  purposes  were  to  achieve  their 
maximum  usefulness  and  employment.1 
For  the  past  decade,  the  U.S.  Atomic 
Energy  Commission  has  sponsored  a 
continuing  research  program  aimed  at 
predicting  amounts  and  extent  of  so-called 
"seismic  damage"  from  underground 
nuclear  detonations.  As  a  result,  theories 


and  methods  have  been  developed  that 
enable  us  to  estimate,  with  reasonable 
accuracies,  the  ground  motion  from 
single  detonations  in  the  sandstone-shale 
medium  most  often  encountered  in  deep 
nuclear  explosions  for  gas  stimulation. 
Ground  motion  parameters  which  appear 
to  best  correlate  with  damage  have  been 
identified,  and  permit  the  amount  of 
damage  to  be  estimated.  Clearly,  a 
number  of  improvements  in  both  motion 
and  damage  predictions  need  to  be  made. 

In  addition,  more  sophisticated  approaches 
have  been  made  which  approach  damage 


1 


predictions  from  a  probabalistic  stand-  present  the  status  of  these  methods.  No 

point.2  In  this  paper  I  attempt  to  review  claim  to  an  exhaustive  treatment  is  made, 

the  steps  by  which  some  of  the  methods  and  it  should  be  kept  in  mind  that  not  all 

now  in  use  have  been  developed,  and  answers  are  presently  in  hand. 


Correlation  of  Damage  with  Peak  Surface  Velocities 


As  an  initial  approach  to  studying 
predictive  methods  of  ground  motion,  we 
investigated  the  experience  of  industry 
with  high  explosives  used  in  quarry  and 

o 

mine  blasting.  We  recognized  from  the 
beginning  that  the  time-history  of  the 
ground  motion  produced  by  a  nuclear 
explosion  would  not  be  the  same  as  that 
produced  from  high  explosives  and  that, 
in  particular,  more  low-frequency 
amplitudes  might  be  expected  from  the 


nuclear  explosion.  However,  the  existing 
blasting  data  seem  to  indicate  that  the 
threshold  for  damage  to  residential 
buildings  (manifesting  itself  primarily  in 
the  cracking  of  plaster)  might  be  inde- 

4 

pendent  of  frequency  of  ground  motion. 
The  data  leading  to  this  conclusion  are 
shown  in  Fig.  1,  which  suggests  that 
below  an  amplitude  of  about  10  cm/sec, 
no  (or  very  little)  damage  might  be  ex¬ 
pected.  Such  a  frequency-independent 


F  =  frequency  in  Hz. 
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damage  behavior  can  also  be  expressed 
in  terms  of  an  energy  ratio  defined  as 
A  /F  ,  where  A  is  the  acceleration  in 
ft /sec  and  F  is  the  frequency  in  Hz. 
Regulations  adopted  by  several  states  of 
the  USA  as  well  as  by  the  U.  S.  Corps  of 
Engineers  set  an  energy  ratio  of  1  as  the 
criterion  for  motion  below  which  damage 
would  not  be  expected.  Such  an  energy 
ratio  of  1  corresponds  to  a  peak  surface 
velocity  of  about  4.8  cm/sec,  assuming 
simple  harmonic  motion. 

This  concept  of  a  damage  threshold 
of  between  5  and  10  cm/sec  was  rudely 
shattered  as  a  result  of  the  Salmon 
explosion,  a  5-kt  detonation  at  a  depth 
of  825  meters  in  a  salt  dome  near 
Hatt i esb ur  g,  Mis s  is s ippi .  App r ox imately 
1000  claims  for  damage  to  buildings  were 
filed  from  the  surrounding  communities.6 
Although  the  few  ground  motion  measure¬ 
ments  made  were  ambiguous,  it  was 
certain  that  surface  velocities  were  very 
much  lower  than  any  damage  threshold 
previously  proposed.  Table  1  lists  Sur¬ 
face  velocities  and  the  percentage  of 
damage  claims  from  three  of  the  com¬ 
munities  in  the  vicinity  of  the  Salmon 
experiment.  Figure  2  shows  an  example 
of  the  type  of  damage  to  buildings  caused 
from  the  ground  motion.  The  type  of 
damage  was  totally  confined  to  what  has 
been  termed  11  architectural"  in  nature. 


Table  1.  Claims  for  building  damage 
from  the  Salmon  experiment. 


Location 

Peak  surface 
velocity 
(cm/sec) 

Claims  /structure 
(%) 

Lumberton 

~1.3 

25 

Purvis 

~0.8 

40 

Hattiesburg 

~0.4 

3 

as  contrasted  with  "structural";  that  is, 
the  damage  involved  cracking  of  interior 
plaster  or  brittle  stucco  or  masonry  walls, 
but  not  structural  elements.  Many  of  the 
cracks  similar  to  the  one  shown  in  the 
figure  closed  with  the  passage  of  time, 
making  a  valid  determination  of  actual 
damage  very  difficult. 

The  Salmon  experience  gave  renewed 
impetus  to  understand  the  mechanism  of 
building  damage  and  to  find  ground  motion 
parameters  which  could  be  correlated 
with  building  damage.  In  particular,  the 
apparent  violation  of  the  previously 
accepted  surface  velocity  criteria  de¬ 
manded  further  study.  This  was  carried 
out  at  the  AECs  Nevada  Test  Site  in 
1965,  in  conjunction  with  the  then  on¬ 
going  series  of  underground  nuclear 
explosions.  At  the  Base  Support  Camp 
for  these  detonations,  a  series  of  new 
buildings  had  been  erected  during  1964. 
These  buildings  were  primarily  one-story 
structures  constructed  of  concrete  block 
and  used  for  housing,  maintenance  facil¬ 
ities,  and  office  space.  Figure  3  shows 
a  typical  building  of  this  type.  Forty- 
three  of  these  buildings  were  selected 
for  detailed  study  and  were  carefully 
examined  inside  and  out  before  and  after 
each  nuclear  explosion,  as  well  as  at  in¬ 
tervals  between  explosions.  Exterior  and 
interior  cracks  and  fractures  were  marked, 
identified,  and  measured,  and  ground 
motion  measurements  were  made  in  the 
area  covered  by  these  buildings.  The 
results  of  this  study6  are  shown  in 
Fig.  4,  in  which  the  total  cumulative 
number  of  cracks  for  all  43  buildings  is 
shown  as  a  function  of  time  in  the  study. 

It  was  found  that,  on  the  average,  2-1/2 
new  cracks  appeared  every  day,  constituting 
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a  natural  cracking  rate  for  these  particular 
buildings.  The  steps  shown  in  Fig.  4  rep¬ 
resent  additional  fractures  caused  by 
nuclear  explosions  taking  place  that  day. 
These  latter  fractures  were  physically  in- 
distinguishable  from  those  which  occurred 
normally.  All  of  them  were  architectural 


in  nature,  and  most  of  them  of  the  type 
easily  overlooked  in  a  casual  inspection. 
An  example  is  shown  in  Fig.  5.  Maximum 
surface  velocity  in  all  cases  was  less  than 
about  0.3  cm/ sec,  confirming  that  some 
type  of  minimal  damage  can  be  caused  by 
motions  of  very  low  amplitude. 


Damage  Correlation  with  Response 
Spectrum  Amplitudes 


It  became  obvious  that  in  order  to 

make  further  progress,  the  amplification 

of  the  ground  motion  by  the  structures 

must  be  considered.  A  natural  avenue  of 

attack  in  this  direction  is  through  the  con- 

7 

cept  of  the  response  spectrum.  This 
method  has  been  used  to  study  natural 
vibrations  of  buildings  and  is  based  on 
the  simplifying  assumption  that  the  re¬ 


sponse  of  a  structure  can  be  approximated 

by  the  envelope  formed  by  the  maximum 

response  of  a  series  of  simple  harmonic 

oscillators  which  are  tuned  to  different 

frequencies  and  damped  to  a  fixed  per- 

8 

centage  of  the  critical  value. 

The  generation  of  the  spectrum  is 
depicted  schematically  in  Fig.  6,  where 
the  maximum  relative  displacement  of 
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quantities  to  be  shown  as  a  function  of 
frequency  or  period  on  a  four-way  graphical 
plot  which  is  especially  convenient  for 
purposes  of  correlating  the  various 
amplitudes  with  damage.  A  typical 
spectrum  is  shown  in  Fig.  7,  which  also 
shows  the  approximate  period  intervals  in 
i  which  buildings  of  various  heights  have 

been  found  to  resonate.  The  relation 
between  spectral  amplitudes  and  resonant 
building  frequencies  is  an  important  one: 
for  example,  if  there  are  no  buildings 
above  three  stories  in  the  area  of  interest, 
the  spectral  amplitudes  above  a  period  of 
about  0.2  sec  are  unimportant  as  far  as 
any  damage  prediction  is  concerned. 

For  one-  or  two- story  buildings, 
where  frequencies  in  the  5  to  20  Hz  range 
(periods  of  0.05  to  0.2  sec)  are  important, 

7 

Nadolski  has  shown  that  the  value  of  the 
pseudo -absolute  acceleration  (PSAA) 
seems  to  correlate  better  with  observed 
damage  than  any  other  motion  parameter. 
The  extensive  ground  motion  measure- 

9 

ments  made  on  the  Rulison  experiment 
(43  ±  8  kt  at  a  depth  of  2570  meters  in 
sand  stone/s  hale)  together  with  thorough 
damage  documentation"^  has  permitted 
this  correlation  to  be  made  much  more 
definitively.  This  has  been  done  by 


Rizer*  and,  more  recently,  by  Farhoomand 
12 

and  Scholl.  Table  2  summarizes  ground 
motions,  spectral  amplitudes,  and  damage 
numbers  from  the  Rulison  experiment. 

These,  together  with  some  additional 
points  from  the  Salmon  experiment  and 
experience  at  Las  Vegas,  are  shown  in 
Fig.  8.  This  figure  shows  that  at  PSAA 
levels  of  close  to  1  g,  in  the  0.05  to  0.2  sec 
period  interval,  over  50%  of  one-  to  two- 
story  structures  can  be  expected  to  sustain 
some  damage.  At  about  0.01  g,  damage 
will  be  sporadic  and  no  (or  very  little) 
damage  can  be  expected  at  lower  values. 

The  line  shown  in  Fig.  8  can  serve  as  a 
quick  and  convenient  method  of  estimating 
building  damage  for  one  or  two  story 
structures  in  the  vicinity  of  nuclear 
detonations,  if  the  response  spectra  at 
these  building  locations  can  be  predicted. 

The  results  from  Rulison  show  that 
even  at  spectral  accelerations  approaching 
the  1  g  level,  actual  damage  was  again 
confined  to  architectural  or  ornamental 
features  of  the  buildings  involved.  Table  3 
summarizes  the  type  of  damage  found  in 
the  communities  surrounding  the  Rulison 
experiment.  Over  half  of  the  total  damage 
involved  cracking  of  interior  plaster  par¬ 
titions  or  damage  to  brick  chimneys. 


Table  2.  Rulison  damage  correlation. 


Location 

Distance 

(km) 

Peak  accel. 
(g) 

PSAA 

0.05 — 0.2  sec 
(g) 

Number  of 
buildings 

Number 

damaged 

Percent 

Grand  Valley 

10.6 

0.55 

0.85 

164 

7  6 

46.5 

Collbran 

19 

— 

0.1 

139 

6 

4.3 

Rifle 

20 

0.10 

0.2 

818 

70 

8.5 

De  Beque 

23 

0.10 

0.15 

106 

6 

5.7 

Silt 

30 

0.035 

0.08 

168 

6 

3.6 

Grand  Junction 

65 

0.017 

0.03 

-4000 

3 

0.075 
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Building  damage 


Fig.  8.  Correlation  between  damage  and  PSAA. 


-10- 


,-V'V: 


While  it  might  be  expected  that  the 
average  repair  costs  would  be  less  at 
lower  motion  levels,  Farhoomand  and 
Scholl  were  not  able  to  establish  such  a 


Examples  of  this  damage  are  shown  in 
Figs.  9-11.  In  all  cases,  the  damage  was 
easily  repaired  at  an  average  cost  of  less 
than  $300  per  claim. 
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Example  of  Rulison  damage 


Fig.  10.  Example  of  Rulison  damage. 
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Fig.  11.  Example  of  Rulison  damage. 
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Table  3.  Types  and  numbers  of  Rulison 
damage. 


Chimneys 

143 

Interior  plaster 

148 

Exterior  walls 

69 

Foundation  walls 

66 

Windows 

24 

Fireplaces 

15 

Furnishings 

22 

Wells  and  cisterns 

27 

Other 

43 

Total 

557 

correlation  for  the  Rulison  damage  ex- 
12 

perience  ;  however,  these  authors  have 
shown  that  a  good  correlation  seems  to 
exist  between  the  maximum  PSAA  in  the 
relevant  frequency  interval  and  the  damage 
repair  cost  divided  by  the  value  of  the 
buildings.  Such  a  procedure,  however, 
may  be  difficult  to  apply  universally  and 
may  also  be  different  for  different  damage 
adjustments  depending,  for  example, 
whether  damage  is  repaired  or  its  cost 
reimbursed. 


Ground  Motion  and  Spectrum  Prediction 


Prerequisite  to  the  use  of  the  method 
presented  above  is  the  prediction  of  the 
response  spectrum  at  the  various  popula¬ 
tion  and  building  centers  which  might  be 
affected  from  a  proposed  nuclear  detona¬ 
tion.  The  prediction  of  peak  values  of 
ground  motion  as  well  as  those  of  response 
spectra  has  been  continuously  refined 
with  the  aid  of  an  extensive  measurements 
program,  both  in  Nevada  and  else¬ 
where.^*^  However,  it  was  not  until 

measurements  from  Gasbuggy  could  be 

1 5 

compared  with  predictions  that  it  be¬ 
came  clear  that  the  scaling  factors 
developed  on  the  basis  of  past  experience 
were  not  adequate  for  very  deeply  buried 
explosions.  Figure  12  compares  Gasbuggy 
observed  and  predicted  peak  accelerations; 
the  measured  values  were  clearly  higher 
than  predicted.  As  a  result,  response 
spectra  were  displaced  towards  higher 
frequencies,  leading  to  higher  PSAA 
values.  These  discrepancies,  as  well  as 

the  detailed  analysis  of  data  from  the 
14 

Nevada  Test  Site,  have  led  Mueller  and 


Murphy  to  a  theoretical  approach  to 
16 

scaling  which  shows  that  both  spectral 
and  ground  motion  amplitudes  depend  not 
only  on  yield  and  medium  properties,  but 
on  the  depth  of  burst  as  well.  The  depth  - 
of-burst  dependence  leads  to  higher  PSAA 
and  higher  surface  acceleration  with 
increased  depth  of  burial.  The  develop¬ 
ment  of  this  theory  uses  the  wave  equation 
and  Hooke's  law  to  relate  the  displacement 
with  the  pressure  function  at  the  elastic 
radius.  The  dependence  on  frequency 
enters  naturally  when  Fourier  transforms 
of  these  functions  are  considered.  For 
two  events  detonated  in  the  same  material, 
the  two  spectra  are  related  by 


z  x  (u>)  | 
^3T| 


I  P1 — el  |  rl-el 
|P2-el|  r2-el 
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where 

|Z(u)  |  =  modulus  of  the  Fourier 

transform  of  the  displace¬ 
ment, 

I  P  ,  I  =  modulus  of  the  Fourier 

transform  of  the  pressure 
at  the  elastic  radius  r^, 

.  where  c  is  the  com- 

0  rel 

press ional  velocity,  and 
|8  s  where  X  and  fj,  are 

the  Lame  constants. 


Both  |  P^|  and  r^  can  be  shown  to  depend 
on  yield  and  overburden  pressure.  By 
using  reasonable  assumptions  for  these 
dependencies,  the  above  equation  can  be 
evaluated  for  a  number  of  yields  and 
depths  of  burst.  Such  a  parameter  study 
can  then  be  used  to  evaluate,  as  a  function 
of  frequency,  the  scaling  exponents  in 
relations  of  the  type 


and 


Z.M  /w.\m<"> 

w  j  t  —  I  I  (for  the  same 
^2^)  y  2/  depth  of  burst) 


Zj(w) 


vn(u) 


(for  the  same 
yield) 


Figures  13  and  14  show  examples  of  m(u>) 
and  n(u)  applicable  to  deeply  buried 
detonations  in  sandstone  and  shale,  such 
as  Rulison. 

Since  in  the  high  frequency  limit  the 
PSAA  amplitude  approaches  the  ground 
acceleration,  Mueller  and  Murphy  derive 
the  following  scaling  law  for  peak  surface 
accelerations: 


When  this  relation  is  used  to  predict 

Rulison  ground  motion  based  on  Gasbuggy 

observations,  the  results  shown  in  Fig.  15 
17 

are  obtained.  This  figure  also  shows 
what  the  prediction  would  have  been  if  the 
depth  correction  had  not  been  applied. 

The  comparison  suggests  that  the  predic¬ 
tion  is  improved  by  including  the  depth 
correction  factor.  The  details  of  the 
spectra  are  more  difficult  to  predict  from 
the  simple  scaling  theory,  as  can  be  seen 
when  the  observed  Rulison  spectra  are 
compared  with  the  predicted  ones;  where 
the  two  disagree,  the  observed  spectra 
are  shifted  to  higher  frequencies  than 
were  expected.  This  situation  is  shown 
in  Figs.  16  and  17.  Some  of  the  dis¬ 
crepancies  are  most  likely  due  to  the 
detailed  geology  of  the  individual  re¬ 
cording  stations.  A  refraction  survey 
at  critical  sites  permits  the  derivation  of 

a  frequency-dependent  amplification 
17 

curve  which  can  be  used  to  correct  the 
spectrum  predicted  for  that  location. 

It  is  instructive  to  compare  some  of 
the  response  spectra  recorded  from 
explosions  with  those  of  a  few  represen¬ 
tative  earthquakes.  This  is  done  in 
Fig.  18.  This  figure  shows  that  the 
maximum  PSAA  at  17  5  km  from  the 
825-kt  Greeley  detonation  is  almost  a 
factor  of  100  below  that  recorded  at  11  km 
from  the  46-kt  Rulison  explosion.  At 
periods  in  excess  of  about  1  sec,  however, 
Greeley  amplitudes  are  very  much  larger 
at  17  5  km  than  the  11  km  Rulison  record," 
Compared  to  Greeley,  the  Rulison  spectrum 
at  a  distance  of  179  km  is  shifted  to  higher 
frequency  and,  of  course,  much  lower 
amplitudes.  Of  the  earthquakes  shown  in 
this  figure,  very  little  damage  was  caused 
by  the  1966  Truckee,  California,  quake. 
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Fig.  14.  Frequency  dependence  of  depth  scaling  exponent. 


1 


while  considerable  damage  resulted 

1 8 

from  the  other  three  shown.  For 
these  earthquakes,  amplitudes  in  the 
resonant  frequencies  of  medium  and  tall 
buildings  are,  of  course,  very  high.  The 
relative  displacements  also  furnish  some 
insight  into  the  degree  of  damage  one 
might  expect  from  either  earthquake  or 


explosion -induced  ground  motion. 

Neither  Greeley  nor  Rulison  caused 
relative  displacements  in  excess  of  about 
1  cm.  At  those  levels,  damage  should 
only  be  architectural  in  nature.  As  dis¬ 
placements  increase,  damage  of  a  more 
serious  nature  becomes  increasingly 
likely. 
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Peak  particle  acceleration 


Fig.  15.  Comparison  of  observed  and  predicted  peak  surface  acceleration — Rulison. 
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5%  pseudo  relative  velocity 
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Fig.  16.  Comparison  of  predicted  and  observed  spectrum — Rulison. 
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5%  pseudo  relative  velocity 


Fig.  17.  Comparison  of  predicted  and  observed  spectrum —  Rulison. 
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